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TECHNICAL MEMORANDUM 


MARSHALL ENGINEERING THERMOSPHERE MODEL, VERSION MET-2007 


1. INTRODUCTION 


The region of the Earth’s atmosphere between about 90 and 500 km altitude is known as the 
thermosphere, while the region above about 500 km is known as the exosphere. For space vehicle 
operations, the neutral atmosphere in these regions is significant. Even at its low density, it pro- 
duces torques and drags on vehicles and affects orbital lifetimes. The thermosphere density above 
100 km altitude also modulates the flux of trapped radiation and orbital debris. Atomic oxygen at 
orbital altitudes is important because it can erode and chemically change exposed vehicle surfaces. 


The NASA Marshall Engineering Thermosphere, version 2007 (MET-2007) is a computer 
program that provides atmospheric density and temperature for the altitude range of 90 to 2,500 
km as a function of latitude, longitude, time, solar flux, and geomagnetic index. Like previous 
versions of the MET model, MET-2007 is a product of the Natural Environments Branch, NASA 
Marshall Space Flight Center (MSFC) and is used for defining the space environment’s neutral 
atmospheric density for vehicle design and mission analysis. MET-2007 retains the capability of the 
previous version (MET-2.0) but with several improvements. These improvements include a correc- 
tion to remove the inconsistency between constituent number density and mass density at certain 
altitudes. Also, minor modifications were made to subroutines using time and date inputs to allow 
for the use of a continuous variable instead of discreet values. 


The sections below provide background on the development of the MET models and 
a description of the model algorithms including the recent updates that were made. 


2. BACKGROUND 


2.1 Jacchia Models 


In the mid-1960s, personnel at MSFC responsible for predicting the orbital decay of satel- 
lites began studies to determine which atmospheric model when combined with the appropriate 
orbit propagation computer programs would most accurately predict the observed altitude decay 
of satellites. Since the observed altitude decay histories were already available, a number of avail- 
able models of the thermosphere were selected for use in appropriate computer programs. Because 
the satellites had decayed prior to the start of the study, actual values of solar proxy input param- 
eters required by the models were used. These proxy parameters were representative of actual solar 
conditions that occurred during the decay periods of the satellites. The Jacchia 1964 (J64) thermo- 
sphere model! when implemented into the obit propagation programs provided the best perfor- 
mance statistically and, therefore, was selected for use by MSFC. 


Almost coincidental with the completion of MSFC’s study, the Smithsonian Astrophysical 
Observatory (SAO) published the Jacchia 1970 (J70) thermosphere model.” The J64 and J70 mod- 
els were similar in their description of the thermosphere but the J70 model included improvements 
based on more recent observations. The major difference between the two models was that the J64 
model temperature-induced density bulge remained on the equator all year while the J70 model 
bulge followed the latitudinal excursions of the Sun, which was supported by observations. Follow- 
on studies showed that the J70 model performed better than the J64 model when implemented 
in obit propagation programs; thus, the J70 model became the baseline thermosphere model used 
by MSFC. 


2.2 Marshall Space Flight Center/Jacchia 70 Orbital Atmosphere Model 


In 1971, SAO published the Jacchia 1971 (J71) model.? Although this model had several 
new features confirmed by observational data, overall, the model was not as representative of the 
atmosphere as the J70 model. An independent study by personnel at MSFC showed that the J70 
model could be improved by adding the J71 model’s seasonal-latitudinal variations in the lower 
thermosphere density below 170 km and adding seasonal-latitudinal variations in helium (He) 
above 500 km altitude while retaining the orbital decay prediction accuracy of the J70 model. This 
new computerized version of the modified J70 model along with the supporting numerical ephem- 
eris and integration algorithms was referred to as the MSFC/J70 Orbital Atmosphere Model.4 


2.3 Marshall Engineering Thermosphere Model 
In the late 1980s, errors were detected in the output of the MSFC/J70 model and were 


traced to the implementation of the J71 modifications described above. These modifications intro- 
duced step function increases in density requiring fairing routines for their elimination. It was also 


determined that there were additional errors attributed to the numerical integration routine. As 

a result, the computer code was rewritten with higher fidelity routines, made more understandable 
and user friendly, and upgraded from FORTRAN IV to FORTRAN 77. This new version of the 
code was named the Marshall Engineering Thermosphere (MET).>7 


While reviewing the 1988 version of MET, it was determined that the model code would 
not correctly adapt to the change in the date of the approaching new millennium. While per- 
forming this update, it was decided to also update the calculation of the solar position to use the 
more recent standard epoch, J2000, which corresponds to January 1.5, 2000, or Julian date (JD) 
2451545.0. An improvement to the model was also achieved in calculating the solar coordinates. 
A low precision ephemeris for the Sun can be calculated using relatively simple equations. This sim- 
pler scheme for computing the solar position requires fewer arithmetic calculations, yet it maintains 
more than adequate accuracy. This modified version of the MET was referred to as MET-1999.8 
This version was also incorporated into the MSFC Global Reference Atmosphere Model (GRAM) 
1999 version. 


A later modification to the MET-1999 version of the model was introduced in the fairing 
subroutine which controls the transition from the lower altitude density profile to that of the upper 
altitude region where the seasonal-latitudinal variation of He is significant. The previous version 
of this subroutine made the transition in a series of small steps. The new scheme introduced in 
MET, version 2.0 (MET-V2.0)? provides a smooth functional form for the transition. 


The latest version of the MET model (MET-2007) incorporates changes discussed in section 4. 
In naming this latest version, it was decided to revert back to the original naming convention of previous 
versions. Since these latest updates were first implemented in the version of the MET model contained 
in GRAM, version 2007, the standalone updated version of MET was named MET-2007. 


3. MODEL DESCRIPTION 


3.1 Overview 


The major source of heating in the thermosphere is from the absorption of the solar 
extreme ultraviolet (EUV) radiation. The temperature in the lower thermosphere increases rapidly 
with increasing altitude beginning at approximately 90 km. Eventually it becomes altitude inde- 
pendent at upper thermospheric altitudes. At these upper altitudes, this asymptotic temperature, 
known as the exospheric temperature, becomes constant with altitude due to the extremely short 
thermal conduction time. The EUV radiation heats only the dayside thermosphere, and although 
conductive and convective processes act to redistribute some of this energy, a large temperature 
gradient always exists between the daytime and nighttime thermosphere. An average daytime 
exospheric temperature is about 1060 K, and an average nighttime exospheric temperature is about 
840 K. 


An additional heat source for the thermosphere is the interaction of the solar wind with 
Earth’s magnetic field. This interaction causes energetic particles to penetrate down into the lower 
thermosphere at high geographic latitudes and directly heat the thermosphere by Joule heating and 
particle precipitation. These energetic particles are also responsible for the aurora seen at these high 
latitudes. 


Whenever the neutral thermospheric gas is heated, it expands radially outward. Because the 
undisturbed thermospheric density decreases with increasing altitude, an outward expansion of 
the gas results in an increase of density at high altitudes. Thus, the daytime thermospheric density 
is greater than the nighttime density, while during times of geomagnetic storms, the high-latitude 
density is greater than it is during undisturbed periods. This anisotropic heating leads to the so- 
called diurnal and polar bulges, which were first inferred from the increased drag experienced by 
orbiting satellites. 


Below the turbopause (located at about 105 km altitude), the atmosphere is well mixed 
by turbulence so that the composition of the atmosphere does not vary with altitude. Above the 
turbopause, however, diffusion becomes so rapid that the altitude variation of the various species 
becomes dependent on molecular weight, with the result that composition varies with altitude. 
Thus, the number densities of the heavier thermospheric species (molecular nitrogen (N5) and 
molecular oxygen (O,)) decrease with increasing altitude much faster than those of the lighter 
species (hydrogen (H) and He). This means that the heavier molecular species predominate in the 
lower thermosphere, while the lighter atomic species predominate in the upper thermosphere. Lift- 
ing of the thermosphere will cause the mean molecular weight at a given altitude to increase, while 
a sinking motion will cause it to decrease. 


As stated above, the NASA MET-2007 model is based upon the J70 model with the inclu- 
sion of equations from the J71 model that characterize the effects of seasonal-latitudinal density 


variations of density below 170 km altitude and that of He above 500 km. Data from satellite drag 
calculations were used to formulate the empirical functions used to calculate the thermosphere tem- 
perature and density. With the proper input parameters, the model calculates an empirically derived 
temperature profile between 90 and 2,500 km altitude for any location on the Earth. With the 
temperature profile specified, the density can be calculated for any altitude and location. The model 
assumes that mixing in the thermosphere prevails to an altitude of 105 km; thus, density between 
90 and 105 km is calculated by the integration of the barometric equation which 1s a function of 
the temperature profile. For altitudes above 105 km, diffusion among the constituents is assumed 

to take place; thus, density is determined by the integration of the diffusion equation which is also 
a function of the temperature profile. 


3.2 Density Variations Associated With Thermosphere Heating 


As stated above, the temperature profile in the J70 model was not the measured quantity 
from satellite drag data. The modeled thermosphere temperature profile was an empirical func- 
tional fit that, when used in the appropriate density equation (barometric or diffusion equation), 
yielded densities that agreed with a database of satellite drag derived densities. 


The thermosphere temperature as a function of altitude is empirically modelled by requiring 
the temperature at 90 km altitude to be constant and equal to 183 K. Then, from 90 km, the tem- 
perature increases to an inflection point at a fixed altitude equal to 125 km. Above 125 km altitude 
the temperature increases asymptotically to an exospheric temperature (7.,). The inflection tem- 
perature and temperature gradient at 125 km altitude is a function of T., as well as the shape of the 
temperature profile. Empirical derived equations using model inputs are used to calculate 7... The 
heating and diurnal variations in thermospheric temperature are modeled by adjusting 7... These 
variations include: 


(1) Variations with solar EUV radiation. 
(2) The diurnal variation. 

(3) Variations with geomagnetic activity. 
(4) The semiannual variation. 


Thus, the exospheric temperature, 7.., is given by 
T,,=T, +ATg +ATg , (1) 


where 77, is the exospheric temperature for a given location on the globe and is a function of the 
EUV heating and diurnal variation, AT¢ is the variation in exospheric temperature associated 
with geomagnetic activity, and AT; is the variation in exospheric temperature due to the observed 
semiannual density variation. A description of these variations is given below. The equations used 
in MET-2007 to calculate the exospheric temperature and its variations as well as the temperature 
profile as a function of altitude are given in references 1 and 8 and summarized in the appendix. 


3.2.1 Solar Extreme Ultraviolet Activity 


The solar EUV radiation that heats the thermosphere is modeled as having two compo- 
nents. One component is related to the active regions on the solar disk and the other is related 
to the overall clear disk itself. The active region component is associated with high-temperature 
solar regions of highly ionized atoms such as active coronal holes, flares, and eruptions. The active 
region component varies rapidly from day to day and is modulated by the 27-day rotation period 
of the Sun. Thus, density variations are often seen to vary with an approximately 27-day periodic- 
ity associated with active regions on the Sun rotating into view. Since EUV measurements are not 
readily measureable from the ground, a surrogate measurement is used that correlates well with 
solar EUV variations. This index is the 10.7 cm (2,800 MHz) solar radio flux (F\9 7) which has been 
measured routinely at ground observatories since 1947 and is an input to the model. 


The EUV radiation from the clear disk component comes from less ionized atoms and 
varies over the course of the 11-year solar cycle. The effect the EUV heating is much less than the 
active region component. For a given increase in the active region component, the thermosphere 
temperature increases by a factor of 3 for the same increase in the disk component in EUV. To 
characterize the disk component of the EUV, Fj, 7 is also used. It was found that the EUV disk 
component is linearly related to the average or smoothed radio flux over a least six solar rotations 
(162 days). Thus, the associated long-term density variation observed over a solar cycle can be 
described quite well by using the smoothed Fj, 7 as an input to the model. 


3.2.2 Diurnal Variation 


Rotation of the Earth induces a diurnal (24-hour period) variation (diurnal tide) in ther- 
mospheric temperature and density. Due to a lag in response of the thermosphere to the EUV heat 
source, density maximizes around 2 p.m. local solar time at orbital altitudes at a latitude approxi- 
mately equal to the subsolar point. The lag decreases with decreasing altitude. Similarly, minimum 
density occurs between 2 and 3 a.m. local solar time at about the same latitude in the opposite 
hemisphere. In the lowest regions of the thermosphere (120 km and below), where characteristic 
thermal conduction time is on the order of a day or more, the diurnal variation is not a predomi- 
nant effect. 


The diurnal variation in exospheric temperature is modelled in MET by assuming that the 
daytime maximum temperature occurs at a latitude equal to the Sun’s declination while the night- 
time minimum temperature, T¢, is diametrically opposed to the maximum temperature at a latitude 
which is minus the Sun’s declination. The diurnal exospheric temperature, 77, at any location on 
the globe is modeled as a function of T¢, latitude, the hour angle of the Sun, and the Sun’s declina- 
tion. The temperature effect of the EUV heating is modelled in the calculation of T; by 


Te = 383+ 3.32Fig.7 +1.8(Fig7- Foz) » 2) 


where Fy, 7 is the 162-day smoothed (six solar rotations) value of the daily average Fj9 7. 


3.2.3 Geomagnetic Activity 


Interaction of solar wind with the Earth’s magnetosphere (referred to as geomagnetic activ- 
ity) leads to a high-latitude heat and momentum source for the thermospheric gases. Energetic 
particle interactions with the Earth’s magnetosphere dissipates energy in the form of Joule heating 
and particle precipitation into the lower thermosphere. This high-latitude heat source is effective 
during both the day and night. Although an intermittent source of energy for the thermosphere, 
it can at times exceed the global EUV energy absorbed by the thermosphere. In addition, though 
the energy is deposited at high latitudes (>60°), the disturbance effects are transmitted to lower 
latitudes through the actions of winds and waves. However, the disturbance effects at low latitudes 
are significantly smaller than they are at higher latitudes. 


The interaction of the solar wind particles with the Earth’s magnetic field is detected as 
disturbances in the geomagnetic field which can be measured by ground-based magnetometers 
and are referred to as geomagnetic storms. Geomagnetic activity is observed to vary with the solar 
cycle, with more intense geomagnetic activity associated with larger solar cycle peaks. 


The heating in the thermosphere associated with variations in geomagnetic activity is mod- 
elled in the MET model by adjusting the exospheric temperature. The variation is given by two 
relationships, one for the geomagnetic index K, and the other for a, (the linearized representation 
of K,) which allows the user of the model to use either index as the input. These indices are based 
on magnetic fluctuation data reported every 3 hours at 12 stations between geomagnetic latitudes 
48° and 63° selected to provide good longitude coverage. The daily planetary geomagnetic index 
(A) is the average of the eight 3-hour a, values measured during a day. There is a time lag between 
changes in the geomagnetic indices and temperature variations which averages approximately 
6.7 hours. Thus, the index that is 6-8 hours before the time of the desired density to be calculated 
should be used in the model. Section 5 describes how to select input values for varying applications. 


3.2.4 Semiannual Variation 


The semiannual variation is not well understood. The fundamental processes driving the 
variation are probably a combination of EUV and Joule heating, meridional transport and mixing, 
variations in lower thermosphere composition, and turbopause height. The variation is observed 
to be latitudinally independent and is modified by compositional effects. Amplitude of the varia- 
tion 1s height dependent and variable from year to year with a primary minimum in July, a primary 
maximum in October, a secondary minimum in January, and a secondary maximum in April. The 
magnitude and altitude dependence of the semiannual oscillation varies considerably from one 
solar cycle to another. This variation 1s important at orbital altitudes and is modeled in the MET 
model as a temperature variation by adjusting the exospheric temperature. The variation in tem- 
perature is a function of the smoothed Fj, 7 and the day of the year. 


3.3 Model Density Seasonal Variations 
There are two seasonal-latitudinal variations in density that are not accounted for by the 


adjustments to the exospheric temperature described above. These are the seasonal-latitudinal 
variations of the lower thermosphere below 170 km altitude and the seasonal-latitudinal variations 


7 


of He above 500 km altitude. These density adjustments are added to the density calculated from 
the barometric or diffusion equation. 


The variations discussed in the previous sections are modelled by adjusting the exospheric 
temperature and keeping the lower boundary temperature fixed at 183 K at 90 km altitude. In real- 
ity, the density at this altitude does vary. The lower thermosphere seasonal-latitudinal variations 
are driven by the dynamics of the lower atmosphere (mesosphere and below). Amplitude of the 
variation maximizes in the lower thermosphere between about 105 and 120 km and diminishes to 
zero around 200 km. Although the temperature oscillation amplitude is quite large, corresponding 
density oscillation amplitude is small. This variation is not important at orbital altitudes. This vari- 
ation is modelled as a function of latitude, day number, and altitude below 170 km. The describing 
equation was adopted from the J71 model. 


Seasonal variations in He measured by satellite mass spectrometers have indicated a strong 
increase in He above the winter pole. Over a year the He number density has been observed to 
vary by a factor of about 40 at 275 km, 10 at 400 km, and 3 or 4 above 500 km. Formation of the 
winter He bulge is primarily due to effects of global scale winds that blow from the summer to 
the winter hemisphere. Amplitude of the bulge decreases with increasing levels of solar activity 
due to increased effectiveness of exospheric transport above 500 km carrying He back to the sum- 
mer hemisphere. This variation is modelled as a function of latitude, Sun declination, and altitude 
above 500 km. The describing equation was also adopted from the J71 model. 


4. MARSHALL ENGINEERING THERMOSPHERE 2007 MODEL UPDATES 


4.1 Corrections to the Marshall Engineering Thermosphere Model Number Densities 
and Mean Molecular Weight 


Output from the MET model includes exospheric temperature, local mass density (p), tem- 
perature (7), pressure (P ), number densities (n,, i= 1,6) for the six MET model constituents (N5, 
O,, atomic oxygen (O), argon (Ar), He, and H), and mean molecular weight (M). If qo(i) denotes 
the individual constituent volume fractions, then number density (n,) is related to mass density by 


n,=qo(i)pN,, 1M, (3) 


where N_,,, is Avogadro’s number. If M; denotes the individual species molecular weights, then the 


mean molecular weight is related to species’ molecular weights and number densities by 


M= > M, ni) Dy (4) 


In examining the MET model constituent output, it was noticed that there were discrep- 
ancies between mass density and number densities as implied by equation (3), or between mean 
molecular weight and number densities as implied by equation (4). These errors are illustrated in 
figure 1, which plots ratios of MET model number density divided by correct number density, or 
MET mean molecular weight divided by correct molecular weight calculated over the globe as a 
function of altitude. Errors in number density of up to 25% are evident below 170 km altitude, 
while errors in mean molecular weight of up to 25% are seen above 440 km altitude. The broaden- 
ing of the curve is due to a dependence of the error on latitude and longitude as well as height. 


Number density discrepancies illustrated in figure 1 were traced to MET subroutine SLY, 
which computes seasonal-latitudinal variation in mass density below 170 km. This subroutine 
adjusted mass density for these effects, but did not adjust number densities to maintain the proper 
relationship of equation (3). This error was corrected in MET by replacing the original number 
densities with new values (1; ) below 170 km, where 


P 6 
eae PNay/ Dp. Mn, ; (5) 


Above 440 km, the error was traced to MET subroutine SLVH, which computes seasonal- 
latitudinal variations in He number density. This subroutine adjusted He number density but did 
not adjust mean molecular weight to retain consistency required by equation (4). This error was 
corrected by replacing the original value of mean molecular weight above 440 km with a new value 
given by equation (4). Note that the mass density from the MET model did not require adjustment 
due to either of these two errors. 


=——— Molecular Weight Ratio 
Number Density Ratio 


Ratio MET/Correct 


102 2 3 4 5 6 7 8 9 103 
Height (km) 


Figure 1. Ratios of uncorrected/corrected molecular weight and ratios 
of uncorrected/corrected number density versus height. 


4.2 Minor Improvements Regarding Model Inputs 


With the implementation of MET-2007 into GRAM-2007, small discontinuities in densities 
at year boundaries, day boundaries, and with time resolutions finer than 1 minute were observed. 
Modifications made to the model to address these issues are listed as follows: 


(1) Treatment of day-of-year as a continuous variable in the semiannual variation term. 
To alleviate discontinuities of density at day boundaries in this routine, the input of a continuous 


variable (day plus day fraction) was implemented. 


(2) Treatment of year as either 365 or 366 days in length (as appropriate), rather than all 
years having a length of 365.2422 days. 


(3) Allowed continuous variation of time input, rather than limiting time increments 
to integer minutes. 
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5. PROGRAM USAGE 


5.1 Input Parameters 


Four basic types of input parameters are required by the model: time, location, solar activ- 
ity, and geomagnetic activity. These may be provided to the model using the sample driver program 
provided, which reads an input data file, puts the data into the INDATA array, and calls the pri- 
mary subroutine MET. The user can choose to replace the driver program or to integrate the MET 
model into a larger program to suit the needs of the project. Note that the INDATA array—and 
therefore all of its elements—is REAL*4. Table 1 presents a summary of the contents of the input 
array. Note that the year is input as a four-digit number, while the month, day, hour, and minute 
are all required to be input as two-digit numbers. The first column contains the array element 
number, the second column identifies the corresponding parameter, and the third column indicates 
the appropriate range of values, where applicable. Note that INDATA 9 controls whether INDATA 


12 is interpreted as a, or K,. 


Table 1. Input parameters. 


Input Variable Parameter 


INDATA (1) | Altitude 


90 to 2,500 km 


Latitude 


-90° to 90° 


Longitude 


-180° to 180° 


Year 


1950 to 2050 


Month 01 to 12 

Day 01 to 31 

Hour 00 to 23 

Minute 00 to 59 

Geomagnetic index type | a,: 2; K,: 1 

Fig 7 0 to 400 x 104 Jansky 
0 to 250 x 104 Jansky 
Geomagnetic index a,: 0 to 400; K,: 0 to 9 


Smoothed F,, 7 


There are three general categories of applications for the MET model: (1) After-the-fact 
calculations of densities, (2) real-time (or near real-time) calculations of densities, and (3) future 
applications in which density estimates are required for some time in the future. In each case, the 
time and location inputs are selected and entered according to the same rules. When entering the 
time, the year is entered in four-digit form, while the month, day, hour, and minute for which the 
calculation is to be done are each entered in two-digit form with the hour and minute being the 
coordinated universal time (UTC). For the position desired, enter the altitude in kilometers and 
geographic latitude and longitude in decimal degrees of the spacecraft location at the time of 
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application. Note that the longitude is measured from the Greenwich meridian eastward. The selec- 
tion of inputs for the solar and geomagnetic activity parameters depends upon which application 
category applies to the calculation as discussed below. 


5.1.1 After-The-Fact Calculations 


This type of calculation is required for times more than 81 days in the past, such as the 
analysis of data from a space mission or testing model performance using observed orbital parame- 
ters for one or more spacecraft. In this case, solar activity is specified using the previous day’s value 
of the Fj, 7 and the centered (about the day for which the calculation is to be done) average of the 
Fo7 over six solar rotations (162 days), referred to as the 162-day smoothed value. The geomag- 
netic index, either a, or K, is the 3-hour index value from 6 to 7 hours prior to the time of applica- 
tion. 


5.1.2 Real-Time Calculations 


When the application requires that the density calculations be made for the current time, 
or less than 81 days in the past or future, averaging the observed Fj, 7 over six solar rotations, as 
described above, is clearly not possible. In this case, the previous day’s value of the Fj 7 1s still used 
for the daily F), 7 value. For the smoothed value, observed values are used where available, and the 
estimated future smoothed value of the Fo7 (typically a 13-month smoothed value) is used as an 
estimate for the balance of the time period where observed values are not available. These estimated 
future values of Fj, 7 are published by the Natural Environments Branch at MSFC and are avail- 
able electronically via the internet (<http://sail.msfc.nasa.gov/>). The observed a, or K, from 6 to 
7 hours prior to the time of application is still appropriate for this application category except for 
the time period up to 81 days in the future, when the estimated future smoothed value of a, is used. 


5.1.3 Future Calculations 


When the calculations of densities must be done for a date more than 81 days in the future; 
e.g., for estimating orbital lifetimes, the problem becomes more dependent upon estimates of future 
average values for the Fj, 7 and the average daily planetary geomagnetic index, A Pp Neither the 
daily nor the 162-day average 9 7 is known from observations, nor is the 3-hour geomagnetic 
activity. In this case, the future estimate of the smoothed value of F9 7 is used as an estimate of 
both the previous day’s value and of the 162-day average value. Similarly, for the input, a future 
estimate of the daily or smoothed index is used as an estimate of the 3-hour geomagnetic index 
input. 
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For future neutral density calculations, the accuracy of the calculation depends primarily 
upon two factors: 


(1) The ability of the thermospheric model to represent the observed neutral density using 
the observed values of solar radio flux (proxy for solar EUV heating) and geomagnetic activity 
used in the development of the model. 


(2) The accuracy with which future solar radio flux and geomagnetic activity can be 
estimated for use as inputs to the thermospheric model. 


5.2 Output Parameters 


Upon completion of the calculation, the MET model passes the results back to the calling 
routine through the OUTDATA array which is REAL*4. These calculated results include the exo- 
spheric temperature, temperature at the input altitude (z), number densities for each constituent, 
the mean molecular weight, total mass density and its logarithm, and total pressure. All of these 
parameters are expressed in mks units. Table 2 details the output array. The first column identifies 
which element of the output array contains the data for the parameter in the second column. The 
total mass density, temperature, and individual species number densities all have the same phase 
variation in the MET model. 


Table 2. Output parameters. 


Output Variable Parameter 


OUTDATA (1) Exospheric temperature (K) 
OUTDATA (2 
OUTDATA (3 
OUTDATA (4 
OUTDATA (5 
OUTDATA (6 
7 


Temperature (K) 


N, number density (m-9) 


O. number density (m-%) 


O number density (m-°) 


Ar number density (m-°) 


) 
) 
) 
) 
) 
) 


He number density (m-3) 


OUTDATA (8) H number density (m-°) 
OUTDATA (9) Average molecular weight 
OUTDATA (10) | Total mass density (kg/m?) 
OUTDATA(11) | Logg mass density (kg/m®) 
OUTDATA(12) | Total pressure (Pa) 


( 
( 
( 
( 
( 
OUTDATA ( 
( 
( 
( 
( 
( 
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APPENDIX A—NUMERICAL EXPRESSIONS FOR TEMPERATURE 
AND DENSITY CALCULATION 


A.1 Expressions for the Temperature Profile 
A.1.1 Temperature Profile Boundary Conditions 


The temperature as a function of altitude T(z) is defined empirically by requiring the 
temperature (7) at altitude zy=90 km to be constant. From 90 km, the temperature increases to 
an inflection point at a fixed altitude z,.= 125 km, having a temperature 7\. and a gradient of G,. 
Above 125 km altitude, the temperature increases asymptotically to an exospheric temperature, T... 
The temperature 7), is a function of T,, and the gradient G,. is a function of T,; thus, also of T,,, 
These boundary values are defined as 


Ty) = 183K (6) 
and 
T,, =a+bT., +cexp(kT.,) : (7) 
where 
a= 444.3807 
b= 0.02385 
c= —392.8292 
k= -0.0021357. 
The temperature gradient is defined as 
- 8 
G, =0 (8) 
and 
T.-T, 
G,,=1.90-~—2 . (9) 
2 2g 


The exospheric temperature is a function of the EUV heating and diurnal variation (77), 
the variation in geomagnetic activity (AT), and the semiannual seasonal variation (AT;), and 
described as 

T,,=T, +ATg+ATz . (10) 


Each of these contributions are described in sections A.1.1.1 through A.1.1.3. 
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A.1.1.1 Diurnal Variation in Exospheric Temperature. The global diurnal exospheric 
minimum temperature is a function of the EUV heating in the upper atmosphere. Fj9 7 is used 
as a proxy index for the EUV variations: 
To = 383+3.32 Fy 7 +1.8(Fig.7-Fig7) - (11) 


The exospheric temperature due to the diurnal variation for a given latitude @ and longitude 
(expressed in terms of the local hour angle H) on the Earth is given by 


cos” n-sin’” @ 


T, =T,(14R, sin” 0) 14+R cos”(t/2) , 12 
- -( us fi 1+R, sin" @ ( ) () 
where 
1 
n= =|6-5, 
1 
O=-|+5, 
T=H+ B+ psin(H+/) 
and 
R,= 0.31 
m =2.5 
n =3 


6, = declination of the Sun (degrees) 
H = hour angle of the Sun (degrees) 


Bp =-3.7° 
pare 
y =43°. 


A.1.1.2 Geomagnetic Variation in Exospheric Temperature. The adjustment to exospheric 
temperature associated with geomagnetic activity is given by two equations, one as a function of 
K, and the other as a function of Ay depending on which geomagnetic index is used as a model 


input: 


ATg = 28K, +0.03 exp(K,] (13) 
and 


AT; =a,+ 100| 1-exp(-0.08a,,) (14) 
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where K, is the 3-hour logarithmic geomagnetic index and dy is the 3-hour linear geomagnetic 
index. There is a time lag of 6.7 hours between changes in the geomagnetic indices and tempera- 
ture. 


A.1.1.3 Semiannual Variation in Exospheric Temperature. The adjustment to the exospheric 
temperature due to the semiannual variation is given by 


AT, = 2.41+ Fio z[ 0.349 + 0.206 (sin(360° t + 226.5°)) |sin(720° t + 247.6°), (15) 


where 


2.16 
c- fons (as{tesn( M4203) 05 (d=number of days since 


January | and Y=length of the tropical year (365.2422 days). 
A.1.2. Temperature Profile 


The temperature profile from 90 < z < 125 km 1s defined as 


T(2)=T,,+>,_j¢n(2-24) 3 (16) 


where 7’, is given by equation (7) and the coefficients c, are given by: 
C= 1.9(T,-Tp)/(z, ~z9); 


c, =0, 


¢,=-1.7(T,-T)/(z,-29) » 


and 


¢4=-0.8(T, -T))/(2,-29) 


The temperature profile for z > 125 km is given by: 


r(e)=1,+ Atan( S2(2-2,)f1+2(2-2,)']), (17) 
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where 


G,, is given by equation (9), 
Ase iT) 
1 


B=4.5x10-° (for z in km) 
n = 2.5. 


A.2 Density Calculations 


A.2.1 Density Between 90 and 105 km Altitude 


Density p(z) at altitudes z between 90 and 105 km is calculated by integrating the baromet- 


ric equation: 


7 T(z) T() 2 M(z)g(z) 
= p| )=—* a a dz 


where 


M = mean molecular weight (function of altitude z) 
T =local temperature (function of altitude z) 
R = universal gas constant (8.31432 J K-! mole7!). 


The values for the parameters at the lower boundary z=z)=90 km are: 


p(Zo) = 3.46x10-? g/cm3 
( 


T (zy) = 183K 


M(zo)= 28.878 g/mole. 


The value of g is calculated by the expression 


g(z)= 


where gy =acceleration of gravity at sea level, 9.80665 m/s”, and R,,=is the Earth radius, 
6.356766 x 10° m. 


(18) 


(19) 


The profile for the mean molecular weight ( M) is given by the empirical equation 


— 6 nN 
M(z)= >) 9 (2-100) (20) 
where the coefficients c,, are: 


Cy = 28.15204 
C1 = 0.085586 

Cy = 1.2840x104 
c, = —1.0056x 10> 
c4 = —1.0210x 10> 
cs= 1.5044x10* 
C6 = 9.9826x10°8. 


The total number density (77) for altitudes between 90 and 105 km is calculated from the 
total mass density p by 


nr =p N,,/M, (21) 
where N_,,, is Avogadro’s number. 


The number density for the constituents N,, Ar, and He is given by 


n, = 0 Gay" ; (22) 
and for O and O,, respectfully, 
n(O) = 2np c rie ; (23) 
0 
and 
M 
n(O,)= nr gilt) -1, (24) 


where the volume fraction g,(i) of the constituents at z=z)=90 km is given as: 


N, = 0.78110 
O, = 0.20955 
Ar = 0.00934 


He = 1.289x10-°. 


A.2.2 Density Above 105 km Altitude 


The density above 105 km is obtained by integrating the diffusion equation to the altitude 
of interest (z) to obtain the number densities of the individual atmospheric constituents, n,. The 
diffusion equation is given by 


l+a 
T(105 M; 
n(z)=n n (105 > ane exp( Ef Sea ) (25) 


where i denotes the individual constituents (N>, O5, O, Ar, He, or H). The parameter a is the 
thermal diffusion factor which is zero for all species except for He. 


For He, the value of —0.380 is used. The parameter M; is the molecular weight of the 
constituent that is being calculated. 


The number density of H is assumed to be negligible below 500 km. Thus, for H, the dif- 


fusion equation is integrated from a lower boundary of 500 km to the altitude of interest. The H 
number density at 500 km is given by: 


2 
logy My (500) = 73.13 - 39.40 logy T,, +5.5(logigT..) . (26) 


where T., is the exospheric temperature as defined in equation (10). The total mass density is 
obtained from the number densities of the constituents by the following: 


6 M; em) (z) 
(A=), 4-3 (27) 
Ny 
thus, the mean molecular weight is given by 
Ny 
in(z)= LO) (28) 


ye ae 
A.2.3 Seasonal-Latitudinal Density Variations 


The lower thermosphere density adjustment for altitudes below 170 km from the J71 model 
is given by 


ca 
4 


Alogy, p(z) = S(z)* Psin’ @, (29) 


where 
= 2 _ | 20 
S(z)=0.014(z -90)exp| -0.0013(z-90)° | and P=sin| = (d+100)}. 
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As before, dis the number of days since January | and Y is the length of the tropical year 
(365.2422 days). 


The He number density adjustment for altitudes above 500 km from the J71 model is given by 


qi) Se os | |, (30) 
4 216, 4 


where € is the obliquity of the ecliptic, 6, is the Sun declination, and @ is the latitude. 


ds 
Alog), 7(He) = 0.65 ra 
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